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We investigate the photon induced tunneling phenomena in a photonic crystal cavity containing a 
strongly coupled quantum dot and describe how this tunneling can be used to generate photon states 
consisting mainly of a particular Fock state. Additionally, we study experimentally the photon- 
induced tunneling as a function of excitation laser power and frequency and show the signature of 
higher manifolds of the Jaynes-Cummings Hamiltonian in the observed photon-statistics. 



A single optical mode confined inside an optical cav- 
ity behaves like a simple harmonic oscillator, where all 
the energy levels are equally spaced. When this cavity 
mode is strongly coupled to a two-level quantum emitter 
such as a quantum dot (QD), the energy structure of the 
coupled system becomes anharmonic. This anharmonic 
(Jaynes-Cummings) ladder has been recently probed in 
atomic [1] and super-conducting [2] cQED system. In ad- 
dition, nonclassical correlations between photons trans- 
mitted through the cavity result from such anharmonic- 
ity, which in turn leads to fundamental phenomena of 
photon blockade and photon induced tunneling. These 
effects have been recently demonstrated in atomic sys- 
tems [3], as well as solid-state [4]. Moreover, photon 
blockade and photon-induced tunneling can be used for 
applications beyond cQED, including generation of single 
photons on demand [5] for quantum information process- 
ing, high precision sensing and metrology [6], as well as 
quantum simulation of complex many-body systems [7] . 
In this Letter, we explore the utility of the photon in- 
duced tunneling and blockade for non-classical light gen- 
eration and probing of higher order dressed states in the 
solid state cQED system consisting of a single quantum 
dot (QD) coupled to a photonic crystal cavity. First, 
we provide numerical simulation data showing that pho- 
ton induced tunneling can be used to preferentially gen- 
erate specific multi-photon states. Following this, we 
present experimental data demonstrating the transition 
from blockade to tunneling regime in such a system and 
show the signature of higher order dressed states ob- 
served in the measured photon statistics. The probing of 
the ladder of dressed states by photon-correlation mea- 
surement has previously been performed experimentally 
only in an atomic cavity QED system [1], while in solid 
state systems it has been studied theoretically [8] and 
signatures of higher order dressed states were observed 
only using four wave mixing [9]. 



The dynamics of a coupled QD-cavity system, coher- 
ently driven by a laser field (Fig. 1 a), is well described 
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FIG. 1: (color online) (a) Schematic of the coupled QD-cavity 
system driven by a Gaussian pulse (coherent state \a)). The 
transmitted light through the cavity is nonclassical(|V>)) due 
to the nonlinearity provided by the strongly coupled QD- 
cavity system, (b) The anharmonic Jaynes-Cummings ladder 
structure. 



by the Jaynes-Cummings Hamiltonian of the form 

H = A a CT + (7_+A c aW.g(a t (7_+acr + )+£(i)(a+a t ) (1) 

which assumes the rotating wave approximation (RWA) 
and a frame of reference rotating with the frequency of 
the laser field cji. Here, A Q = uj a — uji and A c = uj c — 
uii are respectively the detunings of the laser from the 
QD resonant frequency oj a and from the cavity resonance 
frequency u c , g is the coherent coupling strength between 

the QD and the cavity mode, £ (i) = <J [10] is the 
slowly varying envelope of the coherent driving field with 
power P(t) incident onto the cavity with field decay rate 
ft, and a is the annihilation operator for the cavity mode. 
If the excited and ground states of the QD are denoted by 
|e) and \g) then ct_ = \g)(e\ and <j + = \e)(g\. Two main 
loss mechanisms in this system are the cavity field decay 
rate k = u) c /2Q (Q is the quality factor of the resonator) 
and QD spontaneous emission rate 7. When the coupling 
strength g is greater than ^ and 7, the system is in the 
strong coupling regime [11-13]. In this regime, energy 
eigenstates are grouped in two-level manifolds with eigen- 
cnergies given by nuj c ± g^/n (for oj a = u> c ), where n is 
the number of energy quanta in the coupled QD-cavity 
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system (Fig. 1 b). The eigenstates can be written as: 



Signatures of the photon blockade and tunneling can be 
detected through photon-statistics measurements, such 
as the second-order coherence function at time delay zero 

g^(0) = < Lt Q A°^ • ff^(O) i s l css (greater) than 1 in pho- 
ton blockade (tunneling) regime, signifying presence of 
single (multiple) photons in the light coming out of the 
coupled QD-cavity system. g^(0) can be experimentally 
measured by Hanbury-Brown and Twiss (HBT) setup, 
where coincidences between the photons are detected [4] . 
Another important statistical quantity is n th order dif- 
ferential correlation function C^ n \0) = (a^ n a n ) — {a [ a) n , 
which provides a clearer measure of the probability to 
create n photons at once in the cavity [1]. Second order 
differential correlation function can also be expressed as 
C( 2 )(0) = [g^(0) — l]"-c> where n c = (a^a) is the average 
intra-cavity photon number. Particularly for a weakly 
driven system (n c <C 1), C^ 2 \0) becomes positive only 
when the probability of two-photon state becomes sig- 
nificant compared to that of a single-photon state, while 
a peak in C^(0) indicates maximum probability of a 
two-photon state inside the cavity. As the driving power 
increases, the peak in C' 2 '(0) shifts towards detunings 
corresponding to maximum probability of exciting higher 
manifolds, as described below and in the supplement. 

Although the photon blockade and tunneling phenom- 
ena can be observed under continuous wave (CW) exci- 
tation in a numerical simulation [10], for practical con- 
sideration it is important to analyze the response of the 
cavity-QD system to a pulsed driving field. In particular, 
the ability to measure the photon statistics of the sys- 
tem's output during the actual experiment is determined 
by the time resolution capabilities of the single photon 
counters in the HBT setup, which in practice do not allow 
for </ 2 ) measurement of a CW-drivcn cavity-QD system. 
A common way to overcome this limitation is to drive the 
strongly coupled cavity-QD system with a train of weak, 
coherent pulses of sufficiently narrow bandwidth [4] . We 
use quantum trajectory method [14, 15] to analyze the 
pulsed driving of the coupled QD-cavity system and find 
the resulting photon statistics [5] . We also investigate the 
effect of pure QD dephasing [16] on the photon statistics 
and observe that, even though the actual value of g^(0) 
is affected due to dephasing, the qualitative nature of the 
g^ 2 \0) spectrum remains same [10]. As the non-classical 
state is collected from the cavity, only the collapse opera- 
tor corresponding to the cavity decay (a) is monitored. A 
histogram is calculated based on the photon counts in the 
cavity decay channel, and probability P(n) for having ex- 
actly n photons in the system is found. The driving term 
£{t) in the Hamiltonian described in Eqn. 1 is assumed 
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to be of the form £{t) = £ exp{—^pi), where £ is the 
peak amplitude of the pulse. We set t p = 24 .4 ps (i.e., 
full width at half maximum - FWHM of 34 ps), which 
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FIG. 2: (color online) Numerically calculated photon statis- 
tics at the output of the QD-cavity system driven by Gaussian 
pulses with duration r v ~ 24 ps. The simulation parameters 
are g = 2tt x 40 GHz, k = 2tt x 4 GHz, and £ a = 2tt x 9GHz; 
pure QD dephasing is neglected, (a) P(n), probability of gen- 
erating an n photon state at the cQED system output as a 
function of laser-cavity detuning A c . (b) Second order auto- 
correlation g (0) as a function of A c . The red dashed line 
shows the expected <?' 2 '(0) for a coherent state, (c) Second 
order differential correlation C' 2 ^(0) as a function of A c . (d) 
C^ 2 '(0) as a function of the laser-cavity detuning A c for dif- 
ferent values of the peak laser field £ /2n (in units of GHz). 
We observe that the peak in the C' 2 ' (0) occurs at A c = 0.7g 
for weaker excitation (where the second order manifold is ex- 
cited resonantly via two photons). However, with increasing 
excitation power, the peak positions shifts towards A c = 0, 
due to excitation of higher manifolds. 

satisfies the narrow-band condition and corresponds to 
our experimental parameters. 

Figure 2 shows the behavior of the system with better 
than current state of the art [17] but achievable exper- 
imental parameters (assuming QD dipole moment of 30 
Debye embedded in a linear three holes defect cavity with 
mode volume ~ 0.7(A/n) 3 ) resulting in g = 2ir x 40GHz 
and k = 2ir x 4 GHz . These parameters can be achieved 
by improving the alignment of the QD to the cavity 
field and optimizing the photonic-crystal cavity fabrica- 
tion process to achieve higher quality factor. The re- 
sults in Fig. 2a show that such a cavity-QD system can 
be employed to deterministically generate selected Fock 
states of high purity at the cavity output, where the par- 
ticular Fock state can be selected by adjusting the de- 
tuning of the drive laser from the bare cavity resonance 
(no pure QD dephasing is included in the simulation). 
The detuning values (l.lg, 0.9<? and 0.7<jr) are different 
from what one intuitively expects from a lossless strongly 
coupled QD-cavity system under CW driving (g, g/V2 
and g/V3, corresponding to the excitation of first, sec- 
ond and third order manifold, respectively) because of 
both the losses and the pulsed driving of the system [4] . 
We note that, in presence of pure QD dephasing, P{n) 
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FIG. 3: (color online) (a) The transmission spectrum of a 
strongly coupled QD-cavity system showing two polaritons. 
(b) Second order coherence function at t = 0, g' 2 '(0) as a func- 
tion of the laser detuning from the empty cavity frequency. 
The system is excited with r p = 24 ps Gaussian pulses, with 
80 MHz repetition frequency. The dashed grey (solid black) 
line results from a numerical simulation based on the system's 
experimental parameters and no (with) QD blinking. The av- 
erage laser power for the measurement is P avg — 0.2nW. For 
the simulations we use a QD dephasing rate 7d/2/T = 1 GHz. 



for n photon states decreases [10]. From the probability 
distribution of the different Fock states we can find the 
wave- function of the overall photon state \ip) — ^c n |n) 

n 

| 2 , the second order coherence function 

J2 n(n-l)P(r, 

<r J (o) 



with P(n) = \c„ i 

7 (2)("r^ _ (V'|a t Q t oo|V'> 
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and second order 
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differential correlation function C^(0) = (ipla 1 aa\ip) — 
(ipla^ a\i/)) 2 = n ( n ~ — ( n P( n ) ) , which we 

n \ n / 

can measure experimentally. Figure 2b shows g^ 2 \0) as 
a function of A c , the laser detuning from the empty cav- 
ity. The dashed line indicates the expected g^ 2 \0) for 
a coherent state. Figure 2c shows C^(0) as a func- 
tion of A c . C( 2 )(0) transitions from negative to positive 
value with decreased detuning at A c ~ 0.9g, thanks to 
the excitation of the second manifold in the ladder when 
two photons are simultaneously coupled into the cavity- 
QD system. Fig. 2d shows (0) function of A c 
for different laser excitation powers. We note that, the 
peak position changes depending on the excitation laser 
power and at lower driving power we observe the peak 
at A c ~ 0-7g, where the second order manifold is ex- 
cited via two photons. With increasing power, the higher 
(third and more) manifolds starts being populated, and 
the peak in C^ 2 \0) subsequently shifts to smaller val- 
ues of detuning. In Fig. 2c, the peak in C^ 2 '(0) is at a 
detuning of A c ~ 0.5g. 

We confirm our theoretical predictions by performing 
experiment with InAs QDs coupled to a linear three hole 
defect GaAs photonic crystal cavity. Details of the fabri- 
cation and experimental setup can be found in Ref. [4]. 
We measure laser transmission through the system (using 
a cross-polarized reflectivity setup [4]) and observe anti- 
crossing between the QD and cavity (by changing tem- 
perature) signifying the system is in the strong coupling 
regime. At resonance, the QD and cavity mix to generate 
two polaritons, seen as two Lorentzian peaks in Figure 
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FIG. 4: (color online) (a) Normalized differential correlation 
function C' 2 '(0) as a function of the laser detuning. The 
dashed red line shows the result of a numerical simulation 
based on the system's experimental parameters, (b) </ 2 ^(0) 
in the tunneling regime (A c = 0) as a function of laser power 
Pavg measured in front of the objective lens. The solid line 
shows the result of numerical simulation including the effects 
of QD blinking, while the inset plots the numerically simu- 
lated g^ 2 \0) in the absence of blinking. For the simulations 
we use a QD dephasing rate 7d/2/T = 1 GHz. 



3a. By fitting the spectrum at resonance we estimate the 
system parameters as k/2tt = 27 GHz (corresponding to 
Q w 6,000) and ,9/2tt = 21 GHz. To drive the cavity- 
QD system, we use a mode-locked Ti-Sapphire laser, that 
generates 3 ps pulses at a repetition rate of f rep = 80 
MHz. These 3 ps pulses arc passed through a monochro- 
mator to elongate the pulse in time domain, which results 
in pulses with approximately Gaussian temporal profile 
of 34 ps FWHM, corresponding to t p = 24.44 ps (as in 
our theoretical analysis). We determine the amplitude 



of the coherent driving field using £ a 



4Tl2Q Tp f rep h 



[10], where P aV g is the average optical power of the pulse 
train measured before the objective lens and g ~ 0.03 [4] 
is the coupling efficiency of the incident light into the cav- 
ity including all the optics losses. For our experimental 
parameters, £ a ss 2tt\/ P a vg{nW) x 9. 3 GHz. The second 
order auto-correlation (0) is measured as a function of 
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excitation laser frequency (Figure 3b) to observe transi- 
tion from photon blockade to photon induced tunneling 
regime. Typical histograms obtained for blockade and 
tunneling are shown in the supplementary material. We 
estimate (0) as the ratio of the coincidence counts at 
zero-time delay and non-zero time delay. 

Following this we calculate the second order differ- 
ential correlation function C^ 2 ^(0) for the couplcd-QD 
cavity system as a function of the laser-cavity detuning 
(Fig. 4a). We observe the transition of C^(0) from 
negative to positive values and the onset of a peak at 
A c ~ 0.5(? corresponding to the excitation of the higher 
order dressed states. Simulations with our system pa- 
rameters are shown by the dashed line in Fig. 4 a). 
As explained before, the peak in C' 2 '(0) does not cor- 
respond exactly to the resonant excitation of the sec- 
ond order manifold via two photon process, because of 
the additional excitation of the higher order manifolds. 
All the measurements are performed at 14 K. We note 
that, in the simulation, g^(0) in the tunneling regime is 
much larger than the experimentally measured value as 
a result of QD blinking, which causes the experimentally 
collected data to be a weighted average of transmission 
through an empty cavity and a cavity with strongly cou- 
pled QD; in other words, blinking effectively squashes 
the g( 2) (0) curve towards g^iO) = 1 [4]. We model the 
blinking behavior of the QD by assuming that during a 
unit time interval the QD is active for a fraction r and 
inactive for (1— r) of the time. Thus the g^ 2 \0) measured 
in the experiment will be a statistical mixture of the co- 
herent photon state (when QD is inactive, i.e., QD-cavity 
coupling g = 0) and the correlated photons from the cou- 
pled QD-cavity system [10]. We obtain good fit to our 
experimental data with r — 0.65. The vertical error bars 
in all the figures are computed from the uncertainties in 
the fit of the histogram data sets. The horizontal error 
bars are given by the uncertainty in the measurement of 
the laser wavelength or the laser power. 

Finally, Fig. 4 b shows g^ 2 '(0) as a function of ex- 
citation laser power in the tunneling regime (A c = 0). 
This data is taken with the same cQED system on a 
different day, when the cavity is red-shifted compared 
to the previous measurements. For this particular ex- 
periment, the QD and the cavity are resonant at 26 K. 



This slightly higher temperature might cause more QD 
dephasing, leading to a worse value of 5 ( 2 )(0) (1.12 as op- 
posed to 1.2 from the previous measurement). Overall, 
(/ 2 )(0) decreases with increasing laser power as expected 
from the intuitive picture of QD saturation at high driv- 
ing power and the numerical simulation. This clearly 
shows that the bunching observed in tunneling regime is 
coming solely from the quantum mechanical nature of the 
QD-cavity system, and not from any classical effect. We 
also observe interesting oscillatory behavior in g^(Q) as 
a function of power. An oscillatory behavior is also ob- 
served in the simulation that includes the effects of QD 
blinking. Without any QD blinking, the simulation re- 
sults show a mostly monotonically decreasing g^> (0) with 
increasing laser power (inset of Fig. 4b). 

Finally, we would like to point out that these measure- 
ments have been performed at the lowest P avg = 0.2nW 
that we can reliably use, corresponding to £ a w 27r x 
4. 2 GHz. This roughly corresponds to the red plot in the 
theoretical Fig. 2d, where the peak in C^(0) is near 
~ 0.5 g. This lower power limit is caused by the limited 
mechanical stability of the cryostat and the low overall ef- 
ficiency with which we can couple the cavity photons into 
the single photon counters in our HBT setup. The time 
needed to perform the second-order coherence measure- 
ment increases quadratically with decreasing P aV g and 
for low powers the cavity drifts out of focus before we 
can collect sufficient number of coincidence counts. 

In summary, we analyzed the photon induced tunnel- 
ing regime in a coupled QD-cavity system and proposed 
a scheme to use this system for multi-photon state gener- 
ation. In addition, we experimentally characterized the 
second order coherence function g^ (0) for a coupled QD- 
cavity system as a function of laser-cavity detuning and 
laser power. Using the experimental results of the pho- 
ton statistics measurement, we find the signature of the 
higher order manifolds of the Jaynes-Cummings anhar- 
monic ladder in the second order differential correlation 
function C (2 )(0). 
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